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ABSTRACT 

Context. LS 5039 is a Galactic binary system emitting high and very-high energy gamma rays. The gamma-ray flux is modulated on the orbital 
period and the TeV lightcurve shaped by photon-photon annihilation. The observed very-high energy modulation can be reproduced with a 
simple leptonic model but fails to explain the flux detected by HESS at superior conjunction, where gamma rays are fully absorbed. 
Aims. The contribution from an electron-positron pair cascade could be strong and prevail over the primary flux at superior conjunction. 
The created pairs can be isotropized by the magnetic field, resulting in a three-dimensional cascade. The aim of this article is to investigate 
the gamma-ray radiation from this pair cascade in LS 5039. This additional component could account for HESS observations at superior 
conjunction in the system. 

Methods. A semi-analytical and a Monte Carlo method for computing three-dimensional cascade radiation are presented and applied in the 
context of binaries. The cascade is decomposed into discrete generations of particles where electron-positron pairs are assumed to be confined 
at their site of creation. Both methods give similar results. The Monte Carlo approach remains best suited to calculation of a multi-generation 
cascade. 

Results. Three-dimensional cascade radiation contributes significantly at every orbital phase in the TeV lightcurve, and dominates close to 
superior conjunction. The amplitude of the gamma-ray modulation is correctly reproduced for an inclination of the orbit of x 40°. Primary 
pairs should be injected close to the compact object location, otherwise the shape of the modulation is not explained. In addition, synchrotron 
emission from the cascade in X-rays constrains the ambient magnetic field to below 10 G. 

Conclusions. The radiation from a three-dimensional pair cascade can account for the TeV flux detected by HESS at superior conjunction in 
LS 5039, but the very-high energy spectrum at low fluxes remains difficult to explain in this model. 

Key words, radiation mechanisms: non-thermal - stars: individual: LS 5039 - gamma rays: theory - X-rays: binaries 



^1. Introduction 

X, 

• LS 5039 was first id entified as a high-mass X-ray binary 



by iMotch et al.l (Il997h . This binary system is composed of a 
massive O type star and an unknown compact object, possi 



bly a young r otation-powered pulsar dMartocchia et al J 12005 



Dubus 2006bl ). LS 5039 was detected as a ver y high-energy ( 



100 G eV, VHE) gamma-ray source by HESS ([Aharonian et al . 
2005) modulated on the orbital period ( Aharonian et alJl2006l) . 



In a leptonic scenario, the gamma-ray emission is produced 
by inverse Compton scattering of stellar photons on energetic 
electron-positron pairs injected and accelerated by a rotation- 
powered pulsar (pulsar wind nebula scenario) or in a relativistic 
jet powered by accretion on the compact object (microquasar 
scenario). Most of the VHE modulation is probably caused 
by absorption of gamma rays in th e intense UV stellar radi- 



ation field set by the massive star dBottcher & Dermerl 12005 
Bednarekll2006l: lDubusll2006al) . 



Pairs produced in the system can upscatter a substantial 
fraction of the absorbed energy into a new generation of gamma 
rays and initiate a cascade of pairs. The radiation from the 
full cascade can significantly increase the transparency of the 
source, particularly at orbital phases where the gamma-ray 
opacity is high (xyy » 1). A one-zone leptonic model ap- 
plied to LS 50 39 explains the lig htcurve and the spectral fea- 
tures at VHE (IDubus et al.ll2008l) . and yet, this model cannot 
account for the flux detected by HESS at superior conjunc- 
tion where gamma rays should be fully absorbed. Pair cascad- 
ing was mentioned as a possible solution for this disagreement 
jAharonian et al.ll2006h . 

The development of a cascade of pairs depends on the am- 
bient magnetic field intensity. If the magnetic deviations on 
pair trajectories can be neglected, the cascade grows along the 
line joining the source to the observer The cascade is one- 
dimensional. In this case, the cascade contribution is too strong 
close to superior conjunction in LS 5039. A one-dimensional 
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cascad e can be ruled out by HESS observations ( Cerutti et al 
2009bh (see the model in ISierpowska-Bartosik & TorreslbOOsI 



for an alternative solution). If the magnetic field is strong 
enough to deviate and confine electrons in the system, pairs ra- 
diate in_andirecti£ns and a three-dimensional cascade is initi- 
ated ( Bednarekll997 ). The development of a three-dimensional 
cascade i n LS 5039 is possible and was investigated by 
Bednarekl hOOd l2007h wi th a Monte Carlo method and by 
Bosc h-Ramon et al. '('2008a') with a semi-analytical method. 

^osch-Ramon et al. (2008a) derived the non thermal emis- 
sion produced by the first generation of pairs in gamma-ray 
binaries. In their model, the density of secondary pairs is av- 
eraged over angles describing the mean behavior of the radi- 
ating pairs in the system. Here, we aim to investigate the de- 
tailed angular dependence in the gamma-ray emis sion from 
pairs in the cascade. In the microquasar scenario 



Bednarek 



finds consistent flux at superior conjunction in LS 5039 
if the emission originates farther along the jet (> 10 
whose direction is assumed to be perpendicular to the orbital 
plane, including the synchrotron losses. The role of three- 
dimensional casca de is revisited her e in the pulsar wind nebula 
scenario (Marasc hi & Treves! Il98ll: [bubus 2006b), where the 
VHE emitter is close to the compact object location. The aim 
of this article is to corroborate HESS observations of LS 5039 
and to constrain the ambient magnetic field strength in the sys- 
tem, using a semi-analytical and a Monte Carlo computation 
methods. The Monte Carlo code used in the following was pre- 
viously applied to the system Cygnus X-1 for similar reasons 
(IZdziarskietalJl2009h . 

The paper is divided as follows. Sect. 2 gives the main con- 
ditions to initiate a three-dimensional cascade in LS 5039. The 
semi-analytical approach and the Monte Carlo code for cascad- 
ing calculations are presented in Sect. 3 and the main features 
of a three-dimensional pair cascade in binaries are discussed in 
Sect. 4. Sect. 5 is dedicated to the full calculation of a three- 
dimensional cascade in LS 5039. The effect of the ambient 
magnetic field intensity is also investigated in this part. The 
conclusions of the article are exposed in the last section. 

In the following, we use the term "electrons" to refer indif- 
ferently to electrons and positrons. 

2. The magnetic field for 3D cascade 

The development of the cascade is dictated by the intensity of 
the ambient magnetic field in the binary environment. The main 
conditions for the existen ce of a three-dim ensional cascades 
have been investigated by iBednarekl (Il997h and are reviewed 
here and applied to LS 5039. 

The magnetic field B must be high enough to locally 
isotropize pairs once created. This condition is fulfilled if 
the Larmor radius of the pair Rl is shorter than the inverse 
Compton energy losses length given by Acooi - -PeCjelje, 
where je - 1/(1 -01^1^ is the Lorentz factor of the electron 
and je = dje/dt is the Compton energy losses. This provides a 
lower-limit for the magnetic field. In the Thomson regime, this 
is given by 
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Fig. 1. This map shows the domain (gray sur- 
face, 'ISOTROPIC') where a three-dimensional isotropic 
cascade can be initiated as a function of the ambient magnetic 
field B and the energy of the electron Eg. This calculation is ap- 
plied to LS 5039 at periastron (orbital separation d ^ Q.l AU). 
The upper-limit is bounded by the black solid line labeled 
'Bmax' and the lower-limit by the gray solid line 'Bnin ■ For 
B > Bmax ('QUENCHED'), synchrotron losses dominate and 
the cascade is inhibited. For B < B„„„ ('ANISOTROPIC') the 
cascade is not locally isotropized and depends on the magnetic 
field structure. The isotropic domain is truncated at VHE as 
the pairs escape from the system (below the dashed line). 



writing 73 = rJ10\ T^a = T^/4Q 000 K and R^^io = 
7?*/lO/?0 the temperature and radius of the companion star, 
and do.i - d/0.1 AU the orbital separation. Usin g the approxi- 
mate formula for Compton energy losses (Blum enthal & Gould 
1970h . the same condition in the extreme Klein-Nishina regime 
holds if 



Br 



> 1.6x 10"^ tI 



.^lio4.'[ln(r67'*,4) + 2.46] G. (2) 



B^>2xlO-''r3'2^:,4<Vo.lG, 



(1) 



If the Larmor radius is compared with the Compton mean free 
path given by ~ l/n*o-,c, where is the stellar photon 
density and cr,£, the Compton cross section, the condition on the 
magnetic field is more restrictive. In the Thomson regime, the 
electron loses only a small fraction of its total energy per inter- 
action, hence Acooi > ■^w- In the Klein-Nishina regime, most of 
the electron energy is lost in a single scattering and Acooi ~ ^^ic- 
Because the cascade occurs mostly in the Klein-Nishina regime 
in gamma-ray binaries, both conditions lead approximatively to 
the same lower limit for the ambient magnetic field. 

In addition to this condition, pairs are assumed to be 
isotropized at their creation site for simplicity. Pairs will be 
randomized if the ambient magnetic field is disorganized. 
Isotropization of pairs in the cascade will also occur due to 
pitch angle scattering if the magnetic turbulence timescale is 
smaller than the energy loss timescale {e.g. if it is on the order 
of the Larmor timescale). For lower magnetic field intensity 
('anisotropic' domain in Fig. [Ij, the cascade remains three- 
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dimensional but then pairs cannot be considered as locally 
isotropized. In this case, the tr ajectories of the particles should 
be properly computed as in e.g. lSierpowska & B ednarek"2005| 
For B $ 10 G, the cascade is one-dimensional (Cerutti et aj, 
2009bl) . 



If the magnetic field is too strong, pairs locally isotropize 
but cool down via synchrotron radiation rather than by inverse 
Compton scattering. Most of the energy is then emitted in X- 
rays and soft gamma rays, i.e. below the threshold energy for 
pair production. The cascade is quenched as soon as the first 
generation of pairs is produced. This condition gives an upper- 
limit for the magnetic field. Synchrotron losses are smaller than 
inverse Compton losses < Eic for 



Bt ^ 163 n,R.,mdo\ G, 



in the Thomson regime and for 



Bkn ^ 4.7 r6'r^,4-R*,io«^o.'i [in(r6r*,4) + 2.46]'^' G 



(3) 



(4) 



in the deep Klein-Nishina regime dBlumenthal & Gouldl 19701) . 

It can be noticed that the most relevant upper-limit for the mag- 
netic field strength is given by the Thomson formula in Eq. (O, 
since high-energy particles {E^ ^ 1 GeV) with BfCN < B < Bj 
can cool down and get into the cascade domain. 

Figure [1] shows the complete domain where a three- 
dimensional 'isotropic' cascade can be initiated in LS 5039, 
combining the lower and upper-limit for B. This domain en- 
compasses plausible values for the ambient magnetic field in 
the system. It is worthwhile to note that for very high-energy 
electrons E^ ^ 45 c/o.iBo.i TeV, where Bq.i = S/O.l G, the 
Larmor radius becomes greater than the binary separation in 
LS 5039 (Fig. [U. In this case, the local magnetic confinement 
approximation of particles is not appropriate anymore. This is 
unlikely to happen in LS 5039 if the VHE emission has a lep- 
tonic origin since HESS observations shows an energy cut-off 
for photons at « 10 TeV. 

3. Computing methods 

Contrary to the one-dimensional case, three-dimensional pair 
cascading cannot be explicitly computed. Nevertheless, it is 
possible to decompose the cascade into successive generations 
of particles. Two different approaches are presented below, one 
based on semi-analytical calculations and the other on a Monte 
Carlo code. In both models, the primary source of gamma rays 
is point-like and coincident with the compact object position as 
it is depicted in Fig.|2] The origin and the angular dependence 
of the primary gamma-ray flux are not specified at this stage. 
These methods are general and could be applied to any other 
astrophysical context involving 3D pair cascading. 



3.1. Semi-analytical 

A beam of primary gamma rays propagating in the direction de- 
fined by the spherical angles 9 and 4> (see Fig.|2l), produces at a 
distance r to the primary source the first generation of pairs. In 
the point-hke and mono-energetic star approximation, the den- 
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Fig. 2. In this figure is depicted the geometric quantities useful 
for three-dimensional pair cascading calculation in y-ray bina- 
ries. The primary source is point-like and coincides with the 
compact object location. The system is viewed at an angle ^ 
by a distant observer The absorption of primary gamma rays at 
the distance r in the {9, 0) direction creates a secondary source 
of radiation, viewed at an angle ^' by the observer. 



sity of electrons and positrons injected per unit of time, energy 
and volume (s^'erg^'cm^^) is 



(1) 



dtdE.d'V 



dN' 



(0) 
ph 



dtdet dO. 



-gyy 



de\ 



(5) 



where dN'"^lldtde\dQ.ph is the density of primary gamma 
rays of energy e\, ^yy the anisotropic pair production ker- 



nel (lBonometto&Reeslll97U iBottcher & Schhckeiseii 11997 : 
Cerutti et aDl2009bl) and Tyy{r) the yy-opacity integrated from 
the source to the position r. This new density of pairs is spa- 
tially extended and anisotropic but is symmetric with respect 
to the line joining the star to the primary source. For a fixed 
stellar radiation field and a given steady source of primary 
gamma rays, pair production provides a continuous source of 
fresh electrons injected in the binary system environment. 

Pairs are supposed to be immediately confined and 
isotropized by the local magnetic field at their creation site. 
The binary vicinity is surrounded by a plasma of isotropic pairs 
cooling via synchrotron radiation and inverse Compton scatter- 
ing. For simplicity, electrons are assumed to have enough time 
to radiate before escaping their site of injection and the advec- 
tion of particles by the massi ve star wind is ign ored although 
this can have some impact ( Bosch-Ramon et alj|2008a) . For a 
1 TeV electron, the radiative cooling timescales in LS 5039 
are f,c w 20 s (inverse Compton, at the compact object loca- 
tion) and tsyn ~ 400 s (synchrotron, for B = 1 G). The max- 
imum escaping timescale is given by the advection time of 
pairs by the stellar wind. Taking a wind termin al velocity Voo « 
2400 km s"' for the massive star in LS 5039 dMcSwain et al 
20041) . tesc - djvoo « 6 X 10^ s » f,c and fj,,„. In the case where 
pairs would escape the system at the speed of light, electrons 
have just enough time to radiate by inverse Compton scattering 
(tesc - d/c ~ 50 s ^ tic). This extreme situation is unlikely 
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Fig. 3. Spatial distribution of the escaping (i.e. including the effect of gamma-ray absorption) VHE photon density (ph s cm ) 
emitted by the first generation of electrons (isotropized) in the cascade as observed by a distant observer in LS 5039 at superior 
(left) and inferior {right) conjunction. These maps show the gamma-ray density in logarithmic scale (common for both maps), 
where bright and dark regions coiTespond respectively to high and low density. Each map is a slice of the 3D gamma-ray emission 
distribution in the plane that contains the observer (whose direction is indicated by the white solid line) and both stars, computed 
with the semi-analytical method. The primary source of gamma rays is isotropic and lies at the compact object location (origin). 
White dashed lines delimit the eclipsed regions (for the primary source and the observer) by the massive star (bright uniform disk). 
The massive star is assumed point Uke and mono-energetic in the calculations of radiative processes. Distances are normalized 
to the orbital separation. 



since pairs are confined by the ambient magnetic field but pro- 
vides a lower limit for the escaping timescale in the system. 
Assuming that t^sc » fic and f 5,,,, is a rather good approxima- 
tion in LS 5039 for the high-energy particles. 

The steady-state particle d istribution in erg 'cm^^sr ' is 



dOinzburg & Svrovatskiilll964 



dN 



(1) 



dEA'VdD.e 



\Ee\ Je, 



(1) 



I dN, 

An dtdE'd^ 



dE' 



(6) 



with Ef, - Eic + Esyn the inverse Compton and synchrotron 
losses and the volume encircling the binary. Note that the 
annihilation of pairs is not considered in this calculation since 
this effect would be important only for pairs that are almost 
thermali zed. Triplet pair pr oduction y + e"^ — > e''^ -1- -1- 
(s ee e.;?. ^Mastichiadis"l99l') is ignored too (see the discussion 



Cerutti et al. 2009b, Sect. 2.1). 



The total inverse Compton radiation produced by the first 
generation of pairs observed by a distant observer is given by 



dN) 



(1) 



dtde\ dQ.1 



-ff 



dN 



(1) 



dNr 



dEpd'Vd^e dtde 



■e'''"dE,d^, 



(7) 



where is the stellar photon densi ty in cm ^, dNjc/d tdei the 
anisotropic inverse Compton kernel ( Dubus et al.l2008 ) and Tyy 



the absorption from the secondary source up to the observer. 
Depending on the relative position of the secondary source, 
the massive star and the observer, inverse Compton emission 
is anisotropic though pairs are isotropic. The secondary source 
is seen at an angle i/r' with cos iff' = -e* ■ Cghs (Fig- ID) so that 

cos tff' = cos if/ cos ((/(,. - 9) - sin \jj sin {\jjr - ff) cos (p. (8) 

In the point-like star approximation, this viewing angle i/r' is 
related to the interaction angle 9q between photons and elec- 
trons such as cos tfr' - - cos 60. Similarly to inverse Compton 
scattering, the total synchrotron radiation produced by the first 
generation of pairs is 



dN 



(1) 



dNy' dNsyn 



dE.d'V, 



(9) 



dtde\dD.e J J dE^d'VdQ.g dtde\ 
with dNsyn /dtdei the synchrotron kernel averaged over an 
isot ropic distribution of pitch a ngles to the magnetic field (see 



Blumenthal & Gouldll970l) . 



e.g. 

This semi-analytical method can be extended to an arbi- 
trary number of generations. By replacing the primary density 
of gamma rays in Eq. (|5]l by the new density of created photons 
Eqs. (I7]l-(|9]l, the second generation of pairs and gamma-rays in 
the cascade can be computed, and so on for the next genera- 
tions. 
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Fig. 4. The full cascade radiation (all generations) computed 
with the Monte Carlo code (black solid lines) and the pri- 
mary injected gamma-ray source (isotropic, dotted line) are 
shown for [{/ - 30° and 150°. The Monte Carlo output (solid 
gray lines) is compared with the semi-analytical calculations 
(dashed gray lines) in the one-generation cascade approxima- 
tion. There is no magnetic field but pairs are still assumed to 
be confined and isotropized. The massive star is point like and 
mono-energetic. 

3.2. Monte Carlo 

We also used a Monte-Carlo code to simulate the development 
of the full electromagnetic pair cascade in the radiation field of 
the star In this calculation the path and successive interactions 
of photons and leptons are tracked until they escape the system 
(in practice until they reach a distance about 10 t imes the binary 



separa tion). This code was previously used bv iZdziarski et al 



(l2009h to model the TeV emission of Cygn us X-1. It i s simi 
lar in scope and capabilities to the code of Bednared ( 1997 ). 
The present code was developped completely independently, 
and most of the random number generation techniques used 
for computing photon path and simulating the interactions are 
very different from those used by Bednarek. Perhaps the most 
important difference is that the Compton interactions are sim- 
ulated without any approximation, even in the deep Klein- 
Nishina regime. Also, in order to reduce the computing time 
required to achieve high accuracy at high energies, we use a 
weighting technique which avoids following every particle of 
the cascade down to low energies. T he results of both cod es 
were compared and found compatible (IZdziarski et al ] l2009l) . 



4. Three-dimensional pair cascade radiation 

For illustrative purpose only, the primary source of gamma rays 
is assumed isotropic in this section. This assumption allows a 
better appreciation of the intrinsic anisotropic effects of the pair 
cascade emission in binaries. Primary gamma rays are injected 
with a -2 (photon index) power-law spectrum at the location of 
the compact object. For simplicity, the massive star is assumed 



here point-like and mono-energetic. More realistic assumptions 
(injection of isotropic electrons, black body and finite size com- 
panion star) are considered for the calculation of the 3D cas- 
cade emission in LS 5039 in the next Section (Sect.|5ll. 

4.1. Spatial distribution of gamma rays in the cascade 

Figure |3] shows the spatial distribution of the first generation 
of escaping TeV gamma rays seen by a distant observer (i.e. 
including the effect of gamma-ray absorption) produced by the 
cascade in LS 5039 at both conjunctions (for an inclination of 
the orbit / - 60°). These maps are computed with the semi- 
analytical approach. The massive star is assumed point-like for 
the computation of radiative processes but eclipses are con- 
sidered. No pairs can be created behind the star with respect 
to the primary source of gamma rays. Also, gamma rays pro- 
duced behind the star with respect to the observer are excluded 
from the overall cascade radiation (see black regions in Fig.[3]l. 
Synchrotron radiation is neglected in this part: pairs radiate 
only via inverse Compton scattering. 

The spatial distribution of gamma rays is extended and is 
not rotationally symmetric about the line joining the two stars 
(contrary to pairs) since the observed inverse Compton emis- 
sion depends on the peculiar orientation of the observer with 
respect to the binary system. No gamma rays are emitted along 
the line joining the star to the observer direction (see Fig. [3] 
right panel) because pairs undergo rear-end collisions with 
the stellar photons (e* • Cohs - !)■ This effect is smoothed if 
the finite size of the massive star is considered. The escaping 
gamma-ray density at inferior conjunction is more important 
than at superior conjunction as TeV photons suffer less from 
absorption. 

4.2. One and multi-generation cascade 

The semi-analytical method is ideal to study the first genera- 
tion of particles in the cascade as it provides quick and accurate 
solutions. In principle, this method can be extended to an arbi- 
trary number of generation but the computing time increases 
tremendously. The Monte Carlo approach is well suited to treat 
complex three dimensional radiative transfer problems. With 
this method, the full cascade radiation (including all genera- 
tions) can be computed with a reasonable amount of time but a 
large number of events is required to have enough statistics for 
accurate predictions. 

Figure|4]gives the escaping gamma-ray spectra at both con- 
junctions in LS 5039. The Monte Carlo output is compared 
with the semi-analytical results in the same configuration as 
in Fig. [3] for i// - 30° and 150°. Both approaches give sim- 
ilar results for the first generation of gamma rays. There are 
slight differences mainly due to statistical and binning effect 
in the Monte Carlo result, particularly at = 30° where the 
absorption is high. The contribution from additional genera- 
tions of pairs to the cascade radiation is of major importance 
as it dominates the overall escaping gamma-ray flux where 
the primary photons are fully absorbed. The Monte Carlo ap- 
proach is needed to compute the cascade radiation where ab- 
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Fig. 5. Cascade radiation emitted by the first generation computed with the semi-analytical method in LS 5039 at periastron 
for - 30°, 60°, 90°, 120° and 150°. Left: The escaping gamma-ray spectrum (solid line) is compared to the pure-absorbed 
(dashed gray line) and injected (isotropic, dotted line) spectra. The radiation from the cascade only is shown on the right panel. 
Synchrotron radiation is ignored and the massive star is point hke and mono-energetic. 



sorption is strong i.e. at superior conjunction. In practice, the 
one-generation approximation catches the main features of the 
full three-dimensional pair cascade calculation elsewhere along 
the orbit. 



4.3. Comparison with one-dimensional cascade 

Three-dimensional cascade radiation presents identical spec- 
tral features to the one-dimensional limit (Cerutti et al. 2009b) 
(Fig. |5]l. Below the threshold energy for pair production, i.e. 
e\ < m^c* /leod -cos^o) with eo the stellar photon energy, 
pairs cool down via inverse Compton scattering in the Thomson 
regime and accumulate at lower energy in a ~ -1.5 photon 
index power-law tail. Above, emission and absorption com- 
pete, giving rise to a dip in the spectrum. At higher energies 
(ei 10 TeV), the gamma-ray production in the cascade de- 
clines due to Klein-Nishina effect in inverse Compton scatter- 
ing and pair production becomes less efficient. 

Three-dimensional cascade radiation has a strong angular 
dependence (Fig. |5]l that differs significantly from the one- 
dimensional case. Figure |6] presents the modulation of the 
TeV radiation from a ID and 3D cascade along the orbit in 
LS 5039 (the one-dimensiona l cascade radiation is calculated 
with the method described in Cerutti et al. I l2009bl) . iBednarekl 
(l2006h found a similar modulation for the 3D cascade radi- 
ation. Both contributions are anti-correlated. Contrary to the 
one-dimensional cascade, the three-dimensional cascade radi- 
ation preserves the modulation of the primary absorbed source 
of gamma rays since pairs do not propagate. Peaks and dips 
remain at conjunctions. In both cases, the cascade radiation 
flux prevails at superior conjunction where the primary flux 
is highly absorbed. Note that a small dip in the ID cascade 
radiation appears at superior conjunction because absorption 
slightly dominates over emission. The 3D cascade contributes 



less (by a factor ^ 3) than the ID cascade to the total TeV flux 
at this orbital phase. 



4.4. The effect of the ambient magnetic field 

Synchrotron radiation has a significant impact on the cascade 
spectrum. Figure Q shows the effects of an uniform ambient 
magnetic field on the cascade radiation for Z? = 0, 3 and 10 G. 
The VHE emission is quenched as synchrotron radiation be- 
comes the dominant cooling channel for electrons produced in 
the cascade (r,c > f^y,,). The large contribution of the cascade 
in the TeV band is preserved if the magnetic field does not 
exceed a few Gauss (see Fig. [1]). Synchrotron radiation con- 
tributes to the total flux in the X-ray to soft gamma-ray energy 
band. These photons do not participate to the cascade as their 
energy does not exceed 100 MeV, which is insufficient for pair 
production with the stellar photons. 

Figure|2]compares also the contribution from the first gen- 
eration of gamma rays with the full cascade radiation. For low 
magnetic field (B ^ 5 G), all generations should be considered 
in the calculation. For higher magnetic field (B > 5 G), the first 
generation of gamma rays dominates the total cascade radia- 
tion. Only a few pairs can radiate beyond the threshold energy 
for pair production and the cascade is quenched. 

A non-uniform magnetic field was also investigated for a 
toroidal or dipolar magnetic structure generated by the massive 
star {i.e. with a 7? ' or R^^ dependence). These configurations 
do not give different results compared with the uniform case. 
Most of the cascade radiation is produced close to the primary 
source (see Sect. 14. II ) and depends mostly on the magnetic field 
strength at this location. 



B. Cerutti et al.: Modeling the three-dimensional pair cascade in binaries. Application to LS 5039 



7 



1 .00 

> 




ti 

(/) / 
/ 

1 


£\ 

\ 

\ 


1 1 1 1 1 

/ 

/ 




\. ~. 
\ 

\ 

1 - 


h- 




' ^ 

1 / 


— 


1 / 




\\ 


A 














Li_ 








7\ 






^ 0.10 

CU 
N 




Jl \ 




/i \ 






DLU. 


4 


\ 


/ 

/ 1 


1 \ 
/ \ 






O 

z 


1 
I* 


/ \ 


/ 1 


/ \ 






0.01 
















0.0 



0.5 1.0 1.5 

orbital phase (perlastron — 0) 



2.0 



Fig. 6. Modulation of the TeV flux produced by a three- 
dimensional (Monte Carlo calculation, black so lid line) and 
one-dimensional (semi-analytical calculation see ICerutti et al. 

gray solid line) cascade in LS 5039 as a function of the 
orbital phase (two full orbits). Synchrotron radiation is ignored 
for the computation of 3D cascade radiation. The primary ab- 
sorbed flux (identical injection as in Fig. |5j i.e. isotropic) is 
shown (dashed line) for comparison. Conjunctions are indi- 
cated by vertical dotted lines. Orbital parameters are taken from 
[Casares et al. (2005) for an inclination ; - 60°. The companion 
star is point like and mono-energetic. 

5. Three-dimensional cascades in LS 5039 

The full cascade radiation calculation is applied to LS 5039 
and discussed below. The black body spectrum and the spa- 
tial extension of the massive star are taken into account in this 
part. The primary source of gamma r ays is comput ed here fol- 
lowing the model described in .Dubus et al.l (l2008h where the 
pulsar is assumed to inject energetic electron-positron pairs 
with an isotropic power-law energy distribution at the shock 
front, expected to lie at the vicinity of the compact star. Taking 
Voo = 2400 kms~', M - 10"^ yr"' for the massive star 



wind (iMcSwain et al.ll2004h . and a pulsar spin-down luminos- 
ity Lp — 10^^ erg s~\ both wind momenta balance at a distance 
r shock ~ O.lt/ from the pulsar Pairs generated by the pulsar 
emit via inverse Compton scattering on stellar photons the pri- 
mary gamma-ray photons. Contrary to the previous section, the 
primary gamma-ray source is highly an isotropic. The orbita l 
parameters of the system are taken from Casar es et al.l (|2005). 
New opti cal observation s of LS 5039 have been carried out re- 
cently bv lAragona et al. I (l2009l) where slight corrections to the 
orbital parameters have been reported, but these do not change 
the results below. 



5.1. TeV orbital modulation 

The shape of the TeV light curve can be explained with a one- 
zone leptonic model (IDubus et al.ll2008l) that combines emis- 
sion and absorption. However, it overestimates the amplitude 
of the modulation (by a factor > 50 for / = 60°). The TeV 




10" 10' 
e, (eV) 



Fig. 7. Eff'ect of the ambient magnetic field on the cascade ra- 
diation. The cascade is computed with the same parameters 
(Monte Carlo approach) as used in Fig.|4]for ip = 90° with an 
uniform magnetic field Z? = (top), 3 and 10 G (bottom). The 
full escaping gamma-ray spectra (all generations, black lines) is 
compared with the one-generation approximation (gray lines) 
and the injected isotropic spectra (dotted line). The companion 
star is point like and mono-energetic. 



flux observed by HESS varies by about a factor 6 with a min- 
imum at the o rbital phases = 0.1 -0.2 and a maximum at 
<p = 0.8-0.9 (lAharonian et al.l 120061) . The radiation from a 
three-dimensional cascade of pairs decreases the amplitude of 
the TeV modulation yet conserves the light curve pattern (see 
Sect. 14.3b . The flux remains minimum at superior conjunc- 
tion (0 X 0.06) and maximum just after inferior conjunction 
((p ~ 0.85). 

The amplitude of the modulation in LS 5039 can be repro- 
duced for an inclination of the orbit / = 40° (Fig. [8] top panel), 
assuming a c onstant energy den sity of cooled particles along 
the orbit as in lDubus et al.l (l2008h . This assumption imples that 
the injection of fresh particles depends (roughly) as d^^. The 
ambient magnetic field is ^ 1 G (if uniform) otherwise emis- 
sion up to 10 TeV cannot be sustained. For higher inclination 
(/ ^ 50°), the flux at superior conjunction is too small to explain 
observations. For lower inclination (/ ^ 30°), the amplitude of 
the light curve becomes too small. If the injection rate of the 
uncooled primary pairs is instead kept constant along the orbit 
(Fig. [8] bottom panel), a lower inclination (/ 5 30°) is required 
to reproduce an amplitude consistent with observations. Then, 
the light curve presents a broad peak centered at » 0.5. The 
profile of the modulation is not explained to satisfaction in this 
case. 

The cascade radiation contributes significantly at every or- 
bital phase and dominates the overall gamma-ray flux close to 
superior conjunction (0 < < 0.15), where the primary flux 
is highly absorbed. The residual flux observed at superior con- 
junction is explained by the cascade. The averaged spectra at 
high and very-high energy are not significantly changed com- 
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Fig. 8. Theoretical integrated flux above 1 TeV (black solid 
line) in LS 5039 as a function of the orbital phase (two or- 
bits) with an inclination of the orbit / - 40° in both panels. 
The cascade radiation contribution (gray solid line) is com- 
puted with the Monte Carlo approach for a constant injection 
of energy in cooled particles (top) and for a constant injection 
of pairs {bottom) along the orbit. The black-body spectrum and 
the finite size of the companion star are taken into account. 
The ambient magnetic field is small (B < 1 G). Theoretical 
lightcurves are binned in phase interval of width A(f> - 0.1 in 
order to compare with HESS observations (data points) taken 



from Aharonian et al. I (l2006 l). Both conjunctions ('Sup.' and 



'Inf.') are indicated with dotted lines. 




Fig. 9. Theoretical gamma-ray spectum in LS 5039 for 'SUPC 
{i.e. averaged over 0.45 < (p < 0.9, gray dashed line) and 
'INFC (0 < 0.45 or > 0.9 , gray solid line) states as de- 
fined in Aharonian et al. I (l2006h and orbit averaged spectrum 
(black solid line) . Comparison with Fermi (b lack data points, 
Abdo et alj l2009h and HESS (red bowties, lAharonian et al " 
2006h observations. 
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Fig. 10. Same as in Fig.[8](to/7 panel) for ; - 60° with a primary 
source of gamma rays above the compact object and perpendic- 
ular to the orbital plane for an altitude z = 2 



pared with the cas e without cascade (Fig.|9l see also Fig. 6 in 
Dubus et al ]l2008h . It should be noted that the ratio between the 
GeV and the TeV flux decreases if a three-dimensional pair cas- 
cading is considered. The cascade contributes more at TeV than 
at GeV energies with respect to the primary source. If spectra 
are fitted with HESS observations, then the flux expected at 
GeV energies is too low to explain observations. In addition, 
this model cannot ac count for the ene rgy cutoff observed by 
Fermi at a few GeV (lAbdo et al h009'). Electrons radiating at 
GeV and TeV energies may have two different origins. An ex- 
tra component, possibly from the pulsar itself (magnetospheric 



or free pulsar wind emission, see ICerutti et al.l l2009at) might 
dominate at GeV energies. 

5.2. Constraint on the location of tine VHE emitter 

The primary gamma-ray emitter position might not coincide 
with the compact object location. One possibility is to imag- 
ine that particles radiate VHE farther in the orbital plane, for 
instance backward in a shocked pulsar wind collimated by the 
massive star wind. In this case, the primary source is less ab- 
sorbed along the orbit and more power into particles is re- 
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Fig. 11. Orbit averaged spectrum of the first generation of 
gamma rays in LS 5039 with a uniform magnetic field B = 0. 1 , 
1, 5, 10 and 100 G. Compariso n with observations f rom X- 
rays to T eV energies: Suzaku (Takahashi et al. "2OO9I), Fermi 
dAbdo et al 2009.) and HESS (.Aharonian et al.. 2006.) bowties. 



quired to compensate for the decrease of the soft photon density 
from the companion star A consistent amplitude could be ob- 
tained if the primary gamma rays originate from large distances 
10 d), but then the TeV light curve shape is incorrectly re- 
produced as the tendency for the main peak is to shift towards 
superior conjunction. Another possibility is to assume that the 
VHE emitter stands above the orbital plane {e.g. in a jet). This 
situation does not differ significantly from the previous alterna- 
tive. For altitudes z > 2 ^ d, the yy-opacity decreases sig- 
nificantly and the escaping VHE gamma-ray flux increases at 
superior conjunction but the TeV modulation is not reproduced 
as well (Fig. [TOb . Regarding observations, it appears difficult 
with this model to push the gamma-ray emitter at the outer edge 
of the system. The primary source should still lie in the vicinity 
of the compact object {i.e. at distances smaller than the orbital 
separation). 

5.3. Constraint on the ambient magnetic field 

The synchrotron radiation produced by secondary pairs can be 
a dominan t contributor to the overall X -ray luminosity as dis- 
cussed by iBosch-Ramon et al. I (l2008allbl) . Figure [TT| presents 
the orbit-averaged spectrum of the first generation of gamma 
rays in LS 5039 with an inclination / - 40°, using the semi- 
analytical approach for various magnetic field intensity. The 
comparison of the expected flux in the 2-10 keV ba nd with the 
recent Suzaku observations dTakahashi et al. 2009) constrains 
the (uniform) magnetic field strength below 10 G. This result 
is in agreement with the development of a three-dimensional 
cascade (see Sect. 2). The one-generation approximation for 
the cascade is good in this case since for high magnetic field 
{B > 5 G), the contribution from extra-generations can be 
ignored (see Sect. 14.4b . Note that the synchrotron peak en- 
ergy emitted by secondary pairs barely changes with increas- 
ing magnetic field (ei ~ 1 MeV, see Fig. [TTT ). This is due to 



the effect of synchrotron losses on the cooled energy distribu- 
tion of the radiating pairs in the cascade. Synchrotron cooling 
dominates over Compton cooling (?„.„ < f,c) at high energies 
and depletes the most energetic pairs in the steady-state distri- 
bution (see Eq.|6|l. In consequence, the mean energy of cooled 
pairs in the cascade diminishes with increasing magnetic field 
(for a fixed stellar radiation field). The non-trivial combination 
of both effects results in a (almost) constant synchrotron peak 
(the critical energy in synchrotron radiation is proportional to 

6. Conclusion 

Three-dimensional pair cascade can be initiated in gamma-ray 
binaries provided that pairs are confined and isotropized by 
the ambient magnetic field in the system. In LS 5039, a three- 
dimensional pair cascade contributes significantly in the for- 
mation of the VHE radiation at every orbital phase. In partic- 
ular, the cascade radiation prevails over the primary source of 
gamma rays close to superior conjunction {i.e. where the yy- 
opacity is high) and gives a lower flux than the ID cascade at 
this phase. The 3D cascade radiation is modulated differently 
compared with the ID cascade and preserves the modulation 
of the primary absorbed flux because the pairs stay localized. 
In addition, the 3D cascade radiation decreases the amplitude 
of the observed TeV modulation. The amplitude of the HESS 
light curve is correctly reproduced for an inclination of / x 40°. 

The ambient magnetic field in LS 5039 cannot exceed 10 G 
(if uniform) or synchrotron radiation from pairs in the cascade 
would overestimate X-ray observations. This is a reasonable 
constraint as most of massive stars are probably non-magnetic, 
even though strong magnetic fields (> 10 G) hav e been mea- 
sured for a few O stars at their surface (see Donati & Land streetl 
2OO9I for a recent review and references therein). The VHE 



emitter should also remain very close to the compact object 
location, possibly at the collision site between both star winds, 
otherwise the TeV light curve shape is not reproduced although 
this does not rule out complex combinations. 

The model described in this paper is not fully satisfying. 
The spectral shape of VHE gamma rays is still not reproduced 
close to superior conjunction. In addition, the light curve am- 
plitude tends to be overestimated except for low inclinations 
but then the shape is not perfect. It remains difficult to explain 
both the shape and the amplitude of the modulation in LS 5039. 
A possible solution would be to consider a more complex in- 
jection of fresh pairs along the orbit or additional effects such 
as adiabatic losses or advection. A Doppler-boosted emission 
in the primary source can also change the spectru m seen by the 



observer, especially around superior conjunction (IDubus et al 



l2010j). The primary source of gamma rays might be extended, 
VHE photons would come from e.g. the shock front between 
the pulsar wind and the stellar wind or along a relativistic jet. 
The development of an anisotropic 3D cascade is not excluded 
as well. Nevertheless, the calculations show that a three dimen- 
sional pair cascading provides a plausible framework to under- 
stand the TeV modulation in LS 5039. 
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